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Abstract 

The processes pp -» VxV 2 + X -> t x i x hh + X, where V* = W ± or Z and £ { 
denotes a lepton, are calculated to 0(a s ). Total and differential cross sections, 
with acceptance cuts imposed on the final state leptons, are given for the 
Tevatron and LHC center of mass energies. Inclusive and exclusive 0-jet and 1- 
jet cross sections are given. The transverse momenta spectra of the leptons are 
significantly enhanced at high p T by the QCD radiative corrections, especially 
at the LHC energy. Invariant mass and angular distributions are scaled up in 
magnitude by the QCD radiative corrections, but are little changed in shape. 
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I. INTRODUCTION 



The production of massive weak boson pairs {ZZ, W~W + , and W^Z) is an important 
physics topic at high energy hadron colliders. Measurements of these processes are vital 
for testing the Standard Model (SM) and probing beyond it. In particular, these processes 
are important for elucidating the electroweak symmetry breaking mechanism and for testing 
the triple weak boson coupling Futhermore, these processes are also backgrounds to 
new physics signals, for example, new heavy particles such as neutral and charged Higgs 
bosons, techni-mesons, extra gauge bosons, squarks, and gluinos can all decay into weak 
boson pairs. It is therefore important to have precise calculations of hadronic weak boson 
pair production. These calculations should include the leptonic decays of the weak bosons 
since the weak bosons are identified via their leptonic decay products. 

The dominant production mechanism for weak boson pairs in hadronic collisions is via 
the quark-antiquark annihilation subprocess qiq2 — > V{Vi. The cross sections for hadronic 
ZZ, W~W + , and W ± Z production were first calculated in Refs. and 0. Tree level 
calculations of V1V2 + n jets have been given for n < 1 and n < 2 in Refs. 0] and ||, 
respectively. ZZ and W~W + can also be produced via the gluon fusion subprocess gg — > 
VV, which proceeds via a quark box loop and is of order a 2 s The cross section 

for the gluon fusion process is significant at supercollider energies due to the large gluon 
luminosity, but it never dominates the qq annihilation cross section. Weak boson pairs can 
also be produced via the vector boson fusion process in which the incoming quarks radiate 
two vector bosons which subsequently scatter off each other |L(|. This process is mainly 
of interest as a source of Higgs bosons, with the Higgs boson appearing as an s-channel 
resonance. Away from the Higgs boson peak, the weak-boson fusion production rate is only 
a small fraction of the basic q\(j2 — » V1V2 production rate. The QCD radiative corrections 
to hadronic ZZ [IT|,12|, W~W + [I3|,14|, and W^Z |T5|JT6[| production have recently been 
calculated for the case of real weak bosons in the final state. 

In this paper, previous calculations of next-to-leading-order (NLO) ZZ [O, W~W + |L3 
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and W^Z |15j production are extended to include the leptonic decays of the W- and/or Z- 
bosons. Since it is the decay products that are observed in an experiment, the inclusion 
of the leptonic decays in the calculation will make it much more useful for comparing with 
experimental data. Futhermore, cuts can now be applied to the final state leptons, thus 
allowing one to mimic the experimental conditions. The angular distributions of the weak 
boson's decay products are especially important because they are effective spin analyzers 
for the vector bosons |L7|]. They are also of crucial importance when one tries to distin- 



guish the various sources of anomalous couplings in the three-boson vertex [1^]. The NLO 
QCD corrections will have a significant effect on the distributions of the final state leptons, 
especially at the LHC center of mass energy. 

The calculations presented here include the leptonic decays of the weak bosons in the 
narrow width approximation. In this approximation, non-resonance diagrams (e.g., the final 
state ee/i/i, which is formed by the production and decay of ZZ, can also be reached by the 
process pp — ► Z — ► ee with the subsequent radiation of a leptonically decaying Z [Z — > fip] 
from the electron line) are not necessary to maintain gauge invariance. The calculations 
are done by using the Monte Carlo method for NLO calculations [IS] in combination with 



helicity amplitude methods |TJ5|]. The Monte Carlo method for NLO calculations is very 
powerful because any number of observables can be calculated simultaneously by simply 
histogramming the quantities of interest, experimental acceptance cuts can easily be im- 
posed on the calculation, and it is also possible to compute the NLO QCD corrections for 
exclusive channels, e.g., pp — > V\V<z + jet. Helicity amplitude methods make cross section 
calculations tractable for processes involving a large number of tree level diagrams, futher- 
more, the leptonic decays of the weak bosons are trivial to implement at the amplitude 
level. 

The remainder of this paper is organized as follows. The formalism used in the calcu- 
lations is described in Sec. II, numerical results for the Tevatron and LHC center of mass 
energies are given in Sec. Ill, and summary remarks are given in Sec. IV. Technical details 
of the calculation have been relegated to an appendix. 



II. FORMALISM 



Next-to-leading-order calculations of hadronic ZZ, W~W + , and W ± Z production have 
been presented for real weak bosons in the final state. These results are extended in this 
section to include the leptonic decays W — > Iv and Z — > ££ (£ — e, fi). This section begins 
with a brief review of the NLO Monte Carlo formalism used in the calculation and concludes 
with a discussion of the incorporation of the leptonic decays. 

The calculations are done using the narrow width approximation for the leptonically 
decaying weak bosons. This simplifies the calculation greatly for two reasons. First of all, 
it is possible to ignore the contributions from non-resonance Feynman diagrams without 
violating gauge invariance. An example of such a diagram is qq — > Z — > ee followed by 
e — > eZ — > e/i/2, yielding the final state ee/i/i, which is the same final state produced by 
qq — > ZZ — > ee/x/L Secondly, in the narrow width approximation it is particularly easy to 
extend the NLO calculations of real weak boson pairs to include the leptonic decays of the 
W- and Z-bosons. 



A. Monte Carlo Formalism 

The NLO calculations of V1V2 production include contributions from the square of the 
Born diagrams, the interference between the Born diagrams and the virtual one-loop di- 
agrams, and the square of the real emission diagrams. (The Feynman diagrams for the 
case of real weak bosons in the final state can be found in the original references.) The 
calculations have been done using a combination of analytic and Monte Carlo integration 



methods ||18|| . The basic idea is to isolate the soft and collinear singularities associated with 
the real emission subprocesses by partitioning phase space into soft, collinear, and finite 
regions. This is done by introducing theoretical soft and collinear cutoff parameters, 6 S and 



8 C . Using dimensional regularization [20], the soft and collinear singularities are exposed 



as poles in e (the number of space-time dimensions is N = 4 — 2e with e a small number). 



4 



The infrared singularities from the soft and virtual contributions are then explicitly canceled 
while the collinear singularities are factorized and absorbed into the definition of the parton 
distribution functions. The remaining contributions are finite and can be evaluated in four 
dimensions. The Monte Carlo program thus generates n-body (for the Born and virtual 
contributions) and (n + l)-body (for the real emission contributions) final state events. The 
n- and (n+ l)-body contributions both depend on the cutoff parameters S s and S c , however, 
when these contributions are added together to form a suitably inclusive observable, all 
dependence on the cutoff parameters cancels. The numerical results presented in this paper 
are insensitive to variations of the cutoff parameters. 



B. Summary of 0(a s ) V\V 2 production followed by leptonic decays 

The formalism for 0(a s ) hadronic weak boson pair production followed by leptonic 
decays of the weak bosons is summarized here. A detailed discussion about the incorporation 
of leptonic decays into a NLO calculation of a real weak boson pair can be found in Ref. |[21|| , 



where the leptonic decay W — *■ ev was incorporated into a NLO calculation of real 
production. Basically, except for the virtual contribution, all the NLO contributions for real 
V1V2 production have the form 



da NLO ( gi q 2 - V l V 2 ) = da Bo ™( qi q 2 -> V 1 V 2 ) 



1 + 



(1) 



where <7 Born is the lowest order Born contribution, Cf = 4/3 is the quark-gluon vertex color 
factor, and a s is the strong running coupling. Thus the leptonic decays can be incorporated 
by simply making the replacement 

da Bo ™( qi q 2 - VM) — da Bo ™( qi q 2 - V X V 2 - 4W 2 ) (2) 

in the formulas for NLO real V\V 2 production. The leptonic decays are particularly easy to 
incorporate when the calculation is done at the amplitude level; the weak boson polarization 
vectors, e^(k), are simply replaced by the V — > £l decay currents, J^k), in the amplitude. 
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Details of the amplitude level calculations for the Born and real emission subprocesses can 
be found in Ref. @. 

The simple replacement described in the previous paragraph does not hold for the vir- 
tual corrections. Rather than undertake the non-trivial task of recalculating the virtual 
corrections for the case of leptonically decaying weak bosons, we have instead opted to use 
the virtual corrections for real on-shell weak bosons which we subsequently decay ignoring 
spin correlations. Neglecting spin correlations slightly modifies the shapes of the angular 
distributions of the final state leptons, but the total cross sections are not altered as long as 
no angular cuts (e.g., rapidity cuts) are imposed on the final state leptons. For realistic ra- 
pidity cuts, the total cross sections are changed by typically 10% when spin correlations are 
neglected. Since the virtual corrections are small (they are typically less than 10% as large 
as the corresponding Born cross section) and the effects of spin correlations are small, the 
overall result of ignoring spin correlations in the virtual corrections is negligible compared 
to the 20% - 30% uncertainty from the parton distribution functions and the choice of the 
scale Q 2 . (Note that spin correlations are included everywhere in the calculations except in 
the virtual contributions.) 

The results for the NLO calculation of pp — > ViV 2 + X — > l x l\l 2 l 2 + X can now be 
summarized. (The same formalism holds for pp collisions with the obvious replacement 
p —>■ p.) The NLO cross section consists of four- and five-body final state contributions: 

a NLO (pp - V X V 2 + X - tMdi + X) = a^ZyiPP - V X V 2 - hhhh) (3) 

+ a 5 hody (pp -> V X V 2 + X -> l x l x l 2 l 2 + X) . 

The four-body contribution is 

^4 N bo° dy (PP ^V X V 2 ^ W/ 2 ) = a hc + 2 fdv d Xl dx 2 (4) 

91,92 



X 



da NLO 

G qi/p ( Xl , m 2 ) g ?2/p (x 2 , M 2 ) ——( qi q 2 -> v 1 v 2 -> e 1 e 1 e 2 e 2 ) + ( Xl <-> x 2 ) 



dv 

where <r hc is the contribution from the hard collinear remnants (see the Appendix for the 



definition), the sum is over all contributing quark flavors, v is related to the center of mass 
scattering angle 9* by v = |(1 + cos#*), X\ and x 2 are the parton momentum fractions, 
G q / P (x, M 2 ) is a parton distribution function, M 2 is the factorization scale, and 



da 



NLO 



dv 



■(gi92-> ViV 2 ->4W 2 ) 



do - 



Born 



X 



1 + C, 



.ltit 2 t 2 J 



(5) 



a s (/-* 2 
2tt 



.{41n(5 s ) 2 + 31n(A) +4 ln(5 s )ln(A 



+ A FC (9 + jjvr 2 + 3 ln(<5 5 ) - 2 ln(<5 s ) 2 ) } 



da virt 

+ -^—(qiQ2 - V 1 V 2 ) B{V X -> ^A) 5(y 2 



'2*2, 



The hat denotes a parton level cross section, a s (/x 2 ) is the strong running coupling evaluated 
at the renormalization scale /i 2 , 5 S is the soft cutoff parameter, s is the square of the parton 
center of mass energy, and Xpc specifies the factorization convention: Xpc = for the 



universal (Modified Minimal Subtraction MS [22|) convention and Xpc — 1 f° r the physical 



(Deep Inelastic Scattering DIS) convention. The virtual contribution f23 |, da vlvt /dv (q%q 2 
KV2), is multiplied by the V\ — > £\£\ and V 2 — > £ 2 £ 2 branching ratios. 
The five-body contribution is 

^5 bod y (PP -> W 2 + X -> ^W 2 + X) = J2j da{ab -> W 2 c -> 4 W 2 c) 

a,6,c 

G a / P (xi, M 2 ) G b /p(x 2 , M 2 ) + (xi *-> x 2 ) dxi dx 2 , 



(6) 



where the sum is over all partons contributing to the three subprocesses q\q 2 — > V{V 2 g — > 
hhhhg, qi9 — > ViV 2 q 2 — > £i£i£ 2 £ 2 q 2 , and gg 2 — > ViV 2 q~i — > £\£\£ 2 £ 2 q\. The 2^5 subprocess 
is labeled by pi + p 2 — > P3 + P4 + P5 + P6 + P7 and the kinematic invariants s^- and are 
defined by s^- = (p, + Pj) 2 and £^ = (pj — Pj) 2 . The integration over five-body phase 
space and otei cfe 2 is done numerically by standard Monte Carlo techniques. The kinematic 
invariants Sy- and t^- are first tested for soft and collinear singularities. If an invariant for 
a subprocess falls in a soft or collinear region of phase space, the contribution from that 
subprocess is not included in the cross section. The squared matrix elements for the Born 



and real emission subprocesses were evaluated numerically via helicity amplitude methods 
as described in Ref. ||. No attempt has been made to antisymmeterize the amplitudes for 
the case of identical fermions in the final state. 



III. PHENOMENOLOGICAL RESULTS 



The phenomenological implications of NLO QCD corrections to V1Y2 production at the 
Tevatron (pp collisions at ^/s = 1.8 TeV) and the LHC (pp collisions at ^/s = 14 TeV) will 
now be discussed. This section begins with a brief description of the input parameters and 
acceptance cuts. 



A. Input Parameters 

The numerical results presented in this section were obtained using the two-loop ex- 
pression for a s . The QCD scale Aq C d is specified for four flavors of quarks by the choice 
of the parton distribution functions and is adjusted whenever a heavy quark threshold is 
crossed so that a s is a continuous function of Q 2 . The heavy quark masses were taken to 
be m& = 5 GeV and m t = 150 GeV. The SM parameters used in the numerical simulations 
are M z = 91.173 GeV, M w = 80.22 GeV, a{M w ) = 1/128, and sin 2 6 W = 1 - (M w /M z f. 
These values are consistent with recent measurements at LEP, the CERN pp collider, and 
the Tevatron |2^j25|j26|| . The soft and collinear cutoff parameters are fixed to 5 S = 10~ 2 and 
5 C = 10 -3 unless stated otherwise. The parton subprocesses have been summed over u, d, s, 
and c quarks and the Cabibbo mixing angle has been chosen such that cos 2 8c = 0.95. The 
leptonic branching ratios are B(W — > ev) = 0.107 and BiZ — > ee) = 0.034 and the total 
widths of the W- and Z-bosons are Y w = 2.12 GeV and T z = 2.487 GeV. A single scale 
Q 2 = My iV2 , where M Vl v 2 is the invariant mass of the ViV 2 pair, has been used for the 
renormalization scale fi 2 and the factorization scale M 2 . 

In order to get consistent NLO results it is necessary to use parton distribution functions 
which have been fit to next-to-leading order. The numerical results were obtained using the 
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the Martin-Roberts-Stirling (MRS) |27j set SO' distributions with A 4 = 215 MeV. The MRS 
distributions are defined in the universal (MS) scheme and thus the factorization defining 
parameter X FC in Eqs. ([5]) and ( |X2| ) should be Xpc — 0- For convenience, the MRS set SO' 
distributions have also been used for the leading order (LO) calculations. 



B. Cuts 

The cuts imposed in the numerical simulations are motivated by the finite acceptance and 
resolution of the detector. The finite acceptance of the detector is simulated by cuts on the 
four-vectors of the final state particles. These cuts include requirements on the transverse 
momentum, p T , and pseudorapidity, rj = lncot(6 l /2), of the charged leptons and on the 
missing transverse momentum, $ T , associated with the neutrino(s). Charged leptons are also 
required to be separated in azimuthal angle-pseudorapidity space, AR = [(A0) 2 + (A?]) 2 ] 1 / 2 , 
where 9 and <fi are the polar and azimuthal angles relative to the beam; this cut is impelled 
by the finite granularity of the detector. The complete set of cuts can be summarized as 
follows. 

Tevatron LHC 
p T (£) > 20 GeV p T (£) > 25 GeV 

i> T > 20 GeV i> T > 50 GeV 

\ri(£)\ < 2.5 \r}(£)\ < 3.0 

AR(£,£)<0A AR(£,£)<0A 



C. NLO Cross Sections 

The dependence of the total cross section on the soft and collinear cutoff parameters 
is illustrated in Fig. [| which shows the total NLO cross section for pp — > ZZ + X — > 
e~ e + pT fi + + X plotted versus 8 S and S c , for y/s = 14 TeV and the cuts described in Sec. IIIB. 
The n- and n+ 1-body contributions are also plotted for illustration (n = 4 for this process). 



9 



The figure shows that the 4- and 5-body contributions, which separately have no physical 
meaning, vary strongly with 5 S and 8 C , however, the total cross section, which is the sum 
of the 4- and 5-body contributions, is independent of 5 S and 5 C over a wide range of these 
parameters. 

The total LO and NLO cross sections for weak boson pair production, with the cuts 
specified in Sec. IIIB, are given in Table 1 for center of mass energies corresponding to the 
present Tevatron (y/s = 1.8 TeV), an upgraded Tevatron (y/s = 3.5 TeV), and the proposed 
LHC (y/s = 14 TeV). The cross sections in Table 1 have been summed over t — e, /j, and 
both charges of the W in the W ± Z process. 

For definiteness and convenience, differential cross sections will be given for the processes 
pp ^ ZZ + X ^e-e + i2-/2 + +X,pp -> W~W + +X -> e-p e v e e + +X, and pp ^W + Z + 
X — > u e e + fj,~ ji + + X. In practice, the final state leptons would be summed over £ — e, fj, 
and both charges of the W would be summed in the WZ process. Thus the distributions 
given here would be scaled up by the appropriate power of 2. The differential cross sections 
include the cuts described in Sec. IIIB. The figures are arranged in two parts, with parts a) 
and b) being the results for the Tevatron (y/s = 1.8 TeV) and LHC (y/s = 14 TeV) energies, 
respectively. Both NLO and LO results are shown. In some distributions, the NLO 0-jet 
exclusive and the LO 1-jet exclusive cross sections are also given. For these exclusive cross 
sections it is necessary to define a jet. A jet will be defined as a final state quark or gluon 
with 

p T (j)>10GeV and \r)(j)\ < 2.5 (7) 

at the Tevatron, and 

p T (7') > 50 GeV and \r)(j)\ < 3 (8) 

at the LHC. The sum of the NLO 0-jet and the LO 1-jet exclusive cross sections is equal to 
the inclusive NLO cross section. 
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D. ZZ production 

The first process to be considered is pp — > ZZ + X — > e~e + /i~ fi + + X. The invariant 
mass distribution of the four leptons is displayed in Fig. ^|; NLO and LO cross sections 
are shown. The NLO corrections are nearly uniform in the invariant mass at the Tevatron 
energy and increase only slightly with the invariant mass at the LHC energy. 

Figure |3| shows the inclusive differential cross section for the lepton transverse momentum 
(all four leptons have been histogrammed, each with the full event weight). At the Tevatron 
energy, the NLO corrections increase slowly with p T (£), whereas at the LHC energy, the 
corrections increase more rapidly with p T {£). The 0-jet and 1-jet exclusive cross sections are 
also shown. At the Tevatron energy, the 0-jet exclusive cross section is slightly larger than 
the LO cross section, while the 1-jet exclusive cross section is much smaller than the LO 
cross section. At the LHC energy, the 1-jet exclusive cross section is much smaller than the 
LO cross section at small values of p T (£), but becomes comparable to the LO cross section 
at large values of p T {£). The decomposition of the NLO cross section into 0-jet and 1-jet 
components shows that the large NLO corrections at high p T (£) are due to contributions 
from 1-jet real emission subprocesses. The 1-jet exclusive cross section becomes a larger 
fraction of the total NLO cross section at higher energies because the contributions from 
qg initial state processes grow with the center of mass energy due to the increasing gluon 
luminosity. 

The 0-jet and 1-jet exclusive cross sections are of course arbitrary since they depend on 
the jet definition. For example, increasing the p T threshold for a jet will suppress the 1-jet 
and enhance the 0-jet exclusive cross sections (the 0-jet and 1-jet exclusive cross sections 
must sum to the NLO cross section). Nevertheless, Figure |] illustrates that for reasonable jet 
definitions, the 1-jet contribution to the inclusive NLO cross section is small at the Tevatron 
energy, but becomes significant at the LHC energy, especially at high p T (£). 

The inclusive differential cross section for the lepton pseudorapidity, 77 = lncot(0/2) 
where 9 is the polar angle of the lepton with respect to the proton direction in the laboratory 
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frame, is given in Fig. f|. The NLO corrections are largest in the central rapidity region. The 

0- jet and 1-jet exclusive cross sections are also shown. Notice that in the ZZ process, the 

1- jet exclusive cross section is small compared to the NLO cross section. For the W~W + 
and WZ processes, the 1-jet exclusive cross section becomes a larger fraction of the total 
NLO cross section. 

The angular distributions of the leptonic decay products contain information on the 
helicities of the vector bosons. These distributions are simplest in the rest frame of the 
individual vector bosons. For the decay of a polarized Z-boson, Z — > e~e + , the angular 
distributions of the e~ in the Z-boson rest frame are 



dT ^ ft , G F M 



(\ z = 0) = ^4(g 2 y + g 2 A )sm 2 e, (9) 



dcos6 27r v^2 

dT . . , 1 . G F Ml 



dcos9 K 2ttV2 



(\ z = ±l) = ^% (g 2 v + g 2 A )hl + C o S 2 6)±2g v g A co S 6 , (10) 
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where 9 is the angle of the e~ with respect to the longitudinal axis and \ z denotes the 
polarization of the Z-boson; X z = and A z = ±l denote the longitudinal and transverse 
polarizations, respectively. Here Gp is the Fermi coupling constant and g v (g A ) is the 
vector (axial vector) coupling of the Z-boson to fermions. Since the Z-boson coupling to 
charged leptons is almost purely axial vector, transversely polarized Z-bosons produce a 
o(l + cos 2 9) distribution for the e~, while longitudinally polarized Z-bosons yield a sin 2 9 
distribution. Figure [5] shows the polar angle distribution of the e~ in the parent Z-boson 
rest frame, measured with respect to the parent Z-boson direction in the ZZ rest frame, 
i.e., cos9 e - = p e - ■ p z where p e - is the unit-normalized three- momentum of the e~ in the 
parent Z-boson rest frame and p z is the unit-normalized three-momentum of the parent 
Z-boson in the ZZ rest frame. The shapes of the distributions in Fig. |^ indicate that the 
transverse polarizations are dominating the cross section. This is to be expected since the 
qq annihilation process produces Z-boson pairs that are primarily transversely polarized, 
especially at large parton-center-of-mass energies. The sharp drops in the distributions near 
cos6L- = ±1 are due to the kinematic cuts. 
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Figure |6| shows the angular correlation between the decay planes formed by the leptons. 
The angle x between the decay planes is defined by 

(p e - x p e +) ■ x p^) 



cosx 



(11) 



\Pe~ X p e + | X p^+l 

where the momentum vectors are defined in the ZZ rest frame. The NLO and LO curves 
have the same shape at both energies. The shapes of these curves are dominated by the 
effects of the kinematic cuts; without cuts the curves are essentially flat |28 . 



E. W W + production 



Attention now turns to the process pp — > W~W + + X — > e~u e u e e + + X. Since there 
are two invisible neutrinos in the final state, much of the final state kinematic information 
is lost and it is impossible to reconstruct the W~W + invariant mass. The best one can do 



is form the transverse cluster mass [f29] defined by 
M 2 T {c4 T )= \jp 2 cT + ml + \i 



PcT + $1 



(12) 



where c is either a single particle or a cluster of several particles. For the W~W + process the 
cluster is c = e~ + e + . The transverse cluster mass distribution is displayed in Fig. 0. The 
NLO corrections are nearly uniform in M T at both the Tevatron and LHC energies. The 
invariant mass of the charged leptons, M(e~e + ), is also unchanged in shape by the NLO 
corrections. 

Figure |8] shows the inclusive differential cross section for the charged lepton transverse 
momentum (both charged leptons have been histogrammed, each with the full event weight). 
The NLO corrections increase slowly with p T (£) at the Tevatron energy, but at the LHC 
energy, they increase very rapidly with p T (£). Similar behavior is observed in the p T spectra 



of the W^-bosons p^JT^j . The 0-jet and 1-jet exclusive cross sections are also shown. At 
the Tevatron energy, the 0-jet cross section is always larger than the 1-jet cross section. At 
the LHC energy, on the other hand, the 0-jet cross section dominates at small p T (£), while 
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the 1-jet cross section dominates at high p T (£). This behavior is similar to that observed in 
the ZZ process, except now the 1-jet component is a larger fraction of the total NLO cross 
section. 

The missing transverse momentum distribution is presented in Fig. |^. The j> T distribution 
begins to fall rapidly when p" T pa M w . The NLO corrections increase with f T and become 
very large, especially at the LHC energy. The 0-jet and 1-jet cross sections are also shown. 
The 0-jet cross section dominates for f T < M w while the 1-jet cross section dominates 
for j> T > M w . The large NLO corrections at high j> T are due to 1-jet events. At LO the 
VF-bosons are back-to-back in the transverse plane, thus when high p T W-bosons decay, 
the decay product neutrinos will also be nearly back-to-back in the transverse plane. The 
missing transverse momentum, which is the vector sum of the neutrino transverse momenta, 
thus tends to be small due to the acollinear cancellation of the neutrino momenta and the p° T 
distribution falls rapidly. On the other hand, when a W~W + event contains a high p T jet, 
the transverse angle between the decay product neutrinos can easily be acute, thus yielding 
a much larger value of $ T . 



Figure |TD| shows the inclusive pseudorapidity distribution of the charged leptons. The 
NLO corrections are once again largest in the central pseudorapidity region. The 0-jet and 
1-jet exclusive cross sections are also shown. The interesting feature to note is that compared 
to the corresponding distribution for the ZZ process (Fig. f|), the 1-jet cross section for the 
W~W + process is a larger fraction of the total NLO cross section. This trend will continue 
for the WZ process. 



F. WZ production 

The final process to be considered is pp W + Z + X — > v e e + fi~fi + + X. This process 
is of special interest because it is sensitive to the WW Z vertex. The discussion here will 
be limited to the case of standard model couplings at the WW Z vertex. A study of this 
process with NLO corrections and anomalous couplings at the WW Z vertex can be found 
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in Ref. BOI. 



The transverse cluster mass [see Eq. ([T2])], where the cluster is c = e + + fi is shown 



in Fig. [II]. The NLO corrections are nearly uniform in M T at both the Tevatron and LHC 
energies. Since there is only one neutrino in the final state, it is possible to reconstruct the 
WZ invariant mass by requiring the invariant mass of the electron plus neutrino system to be 
equal to the VT-boson mass. This constraint gives a quadratic solution for the longitudinal 
momentum of the neutrino and thus there is a two-fold ambiguity in the reconstructed 
WZ invariant mass. The reconstructed WZ invariant mass, formed by histogramming both 
reconstructed invariant masses, each with half the event weight, is qualitatively similar to 
the cluster transverse mass. In particular, the shape of the reconstructed WZ invariant 
mass is unchanged by the NLO corrections. 

Figure shows the inclusive transverse momentum distribution of the charged leptons 
(all three charged leptons have been histogrammed, each with the full event weight). The 
NLO corrections once again increase with p T (£) and are especially large at the LHC energy. 
The 0-jet exclusive cross section is the dominant component of the inclusive NLO cross 
section at the Tevatron energy, whereas at the LHC, the 0-jet cross section dominates at 
small Pt{£) while the 1-jet cross section dominates at high p T (£). The j) T distribution, which 
is shown in Fig. [L3[ exhibits the same qualitative features as the Pt{£) distribution. 



The inclusive pseudorapidity distribution of the charged leptons is displayed in Fig. [14. 
Note that the 1-jet cross section is now an even larger fraction of the total NLO cross section 
than it was in either the ZZ or W~W + process (see Figs. |] and pT0|) . 



Figure [15| shows the polar angle distribution of the e + in the VT-boson rest frame, mea- 
sured with respect to the H^-boson direction in the W + Z rest frame, i.e., cos6 e + = p e + ■ p w 
where p e + is the unit-normalized three-momentum of the e + in the H^-boson rest frame and 
p w is the unit-normalized three-momentum of the H^-boson in the WZ rest frame. The 
cos 9 e+ distribution, and the analogous cos 9^- distribution discussed in the next paragraph, 
both contain two-fold ambiguities corresponding to the two solutions for the longitudinal 
momentum of the neutrino. Both solutions have been histogrammed, each with half the 

15 



event weight. For the decay of a polarized VT-boson, W + — > e + u e , the angular distributions 
of the e + in the P^-boson rest frame are 

(X w = 0) = — — sin , (13) 



dcos0 v ^ 87ta/2 



(A w = ±l) = ^— ^-(licos^, (14) 
a cos 87TV 2 2 

where # is the angle of the e + with respect to the longitudinal axis and X w = (±1) denotes 
the longitudinal (transverse) polarization(s) of the H /+ -boson. The qq annihilation process 
produces WZ pairs that are primarily transversely polarized, futhermore, for qiq~2 — > W + Z, 
the helicity combination (X w = — 1, X z = 1) gives the dominant contribution pi ]. This 
explains the |(1 — COS6 1 ) 2 shape of the cos9 e + distribution in Fig. [15]. The distributions fall 
near cos# e + = ±1 due to the kinematic cuts. 

Figure shows the polar angle distribution of the yT in the Z-boson rest frame, mea- 
sured with respect to the Z-boson direction in the W + Z rest frame, i.e., cos 6^- = p^- ■ p z 
where p„- is the unit-normalized three-momentum of the [i~ in the Z-boson rest frame and 
p z is the unit-normalized three-momentum of the Z-boson in the WZ rest frame. Vestiges 
of the |(1 + cos 6) distribution characteristic of transversely polarized Z-bosons can be seen 
in the figure. The asymmetry in the distributions comes from the cos# term in Eq. (|1(]) and 
has a negative slope because the dominant contribution comes from the helicity combination 
(X w = -l,X z = l). 



G. Discussion 



Comparing the three processes, one sees that the NLO corrections increase for the pro- 
cesses in the order ZZ , W~W + , WZ. The 1-jet exclusive cross section also becomes a larger 
fraction of the total NLO cross section in this same order. Both of these features are due to 
the qg — > ViV^q real emission subprocesses. 

A decomposition of the NLO corrections into components from order a s qq and qg initial 
states shows that the qq components are of similar size in all three processes (the qq com- 
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ponent is approximately 10% the size of the corresponding Born cross section), whereas the 
order a s qg component increases for the processes in the order ZZ, W~W + , WZ 



A similar decomposition of the 1-jet exclusive cross section exhibits the same behavior. 

The NLO corrections and the 1-jet fraction are largest for the WZ process because the 
Born cross section for this process is suppressed due to destructive interference between 
the S-, T-, and [/-channel Feynman diagrams. This destructive interference produces an 
approximate amplitude zero |3lj ; the dominant helicity amplitudes (Xw = ±l,Az = =Fl) 



have an exact zero, whereas the other helicity amplitudes remain finite but small. The 
situation is analogous to the qq —> Wj process which, because of the massless photon, 
has an exact zero in all the helicity amplitudes fl32| . The NLO QCD corrections to the 
process qq — > W'j are very large at high center of mass energies as a result of the destructive 
interference and the large gluon luminosity f33| . The approximate amplitude zero in the 



qq — > WZ process produces dips in the distributions of cos6** [[H]], y* z , and y z — y w |16 



where 9* is the center of mass scattering angle, y* z is the Z-boson rapidity in the WZ center 
of mass frame, and y z (y w ) is the Z-boson (W-boson) rapidity in the laboratory frame. The 
approximate amplitude zero also suppresses the WZ Born cross section at high p T (Z). The 
qg — > WZq subprocess, on the other hand, is not suppressed by destructive interference. As 
a result, the NLO corrections and the 1-jet fraction are larger for the WZ process than for 
either the ZZ or W~W + process. 

In all three processes, the NLO corrections are largest at high p T because the qg — » V-^Viq 
subprocesses are enhanced in diagrams where a weak boson and a quark are produced at 
high p T , with the quark radiating the other weak boson, i.e., diagrams in which qg — > V\q 
followed by q — > qV^. These subprocesses are enhanced by a factor log 2 (p T (Vi) / My 2 ) []TB| , p~4[] , 
which arises from the kinematic region where Vi is nearly collinear to the quark. 

The W~W + process has a larger NLO correction and a larger 1-jet fraction than the 
ZZ process. These features are again due to differences in the relative importance of the qg 
initial state processes. The qg initial state processes are of more relative importance for the 
W~W + process than for the ZZ process, however, it is not clear why this is so. The mass 
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difference between the W- and Z-boson is not enough to explain the difference. 

Another notable feature of the NLO corrections is that they increase with the center of 
mass energy. This behavior is also due to the qg initial state processes. The contributions 
from these processes increases with the center of mass energy due to the increasing the gluon 
luminosity. 

IV. SUMMARY 

The QCD radiative corrections to hadronic ZZ, W~W + , and W ± Z production have 
been calculated to order a s with leptonic decays of the weak bosons included. The inclusion 
of the leptonic decays makes the calculations more realistic since it is the leptonic decay 
products that are observed in an experiment. Distributions of the final state decay products 
have been given for both inclusive and exclusive channels for the Tevatron and LHC center 
of mass energies. The calculations include typical acceptance cuts on the final state leptons. 

The calculations were done by using the Monte Carlo method for NLO calculations in 
combination with helicity amplitude methods. With the Monte Carlo method it is easy 
to impose experimentally motivated acceptance cuts on the final state leptons, also, it is 
possible to calculate the order a s QCD corrections for exclusive channels, e.g., pp — > ViV 2 + 
jet. The narrow width approximation has been used for the decaying weak bosons. This 
simplifies the calculation greatly since it is possible to ignore contributions from non-resonant 
Feynman diagrams without violating gauge invariance. Futhermore, in the narrow width 
approximation it is particularly easy to extend previous NLO calculations of real weak boson 
pairs to include the leptonic decays of the weak bosons. Spin correlations are included 
everywhere in the calculation except in the virtual contributions where they can be safely 
neglected. 

The QCD radiative corrections enhance the transverse momenta spectra of the leptons 
at high p T) especially at the LHC energy. These enhancements are due to the opening of 
the qg subprocesses at order a s . The contribution from these processes increases with the 
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center of mass energy due to the increasing gluon density in the proton. Invariant mass 
and angular distributions are scaled up in magnitude by the QCD radiative corrections, but 
undergo little change in shape. 



The author would like to thank U. Baur, J. F. Gunion, T. Han, and D. Zeppenfeld for 
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The real emission subprocesses qi(pi) + #2(^2) — > V^Vi9 — > ^1(^3) + h(P4) + ^2(^5) + 
h(pe) + g{P7) (and the cross subprocesses) have hard collinear singularities when tn — > or 
t 2 7 — > [tij = {pi — Pj) 2 ]- These singularities must be factorized and absorbed into the initial 
state parton distribution functions. After the factorization is performed, the contribution 
from the remnants of the hard collinear singularities has the form 
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APPENDIX: HARD COLLINEAR CORRECTIONS 




da 



(Al) 




1-<5 S 



dz 



z 



G, 2 /p(AM 2 )P qq (z) 






with 




(A2) 
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The Altarelli-Parisi splitting functions in A^ = 4 — 2e dimensions for < z < 1 are 

~l + z 2 



1 - Z 



el-* 



P qg (z,e) = —- 



z 2 + (l- z? - e 



(A3) 



(A4) 



and can be written as 



P tj (z,e) = P ij (z) + eP^(z) 



(A5) 



which defines the PL functions. The functions F qq and F qg depend on the choice of factoriza- 
tion convention and the parameter Xpc specifies the factorization convention; Xpc — for 
the universal (Modified Minimal Subtraction MS |p2| ) convention and Xfc = 1 for the physi- 
cal (Deep Inelastic Scattering DIS) convention. For the physical convention the factorization 
functions are 



F qq (z) 



c F 
1 



i + z 2 . n 

In 

1 - z 



{z 2 + (l-z) 2 }ln 



2 1 - z 



+ 2z + 3 



+ 8z(l-z) 



1 



(A6) 
(A7) 



The transformation between the MS and DIS schemes is discussed in Ref. || . The param- 
eter M 2 is the factorization scale which must be specified in the process of factorizing the 
collinear singularity. Basically, it determines how much of the collinear term is absorbed 
into the various parton distribution functions. 
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TABLES 

TABLE I. Total cross sections for pp -► ZZ + X -► ij&h + A", 
-► TT"VL+ + X ^1^2 + X, and -► VLZ + X -► 4M2^2 + A" for center of mass 
energies corresponding to the present Tevatron, an upgraded Tevatron, and the proposed LHC. 
The cross sections have been summed over £ = e, p and both charges of the W in the WZ process. 
The cuts listed in Sec. IIIB have been imposed. 



(TeV) 


(-) 
pp 




a(ZZ) (fb) 


a(W'W + ) (fb) 


a(W ± Z) (fb) 


1.8 


pp 


LO 
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190. 
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NLO 


6.2 


260. 


22. 
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pp 


LO 


12. 


440. 


44. 






NLO 


14. 


590. 


59. 


14. 


pp 


LO 


36. 


570. 


43. 






NLO 


43. 


960. 


77. 
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FIGURES 

FIG. 1. The dependance of the total NLO cross section for pp — > ZZ + X — > e~e + p,~[i + + X 
at ^fs = 14 TeV on the soft and collinear cutoff parameters. In part a) the total NLO cross section 
is plotted versus the soft cutoff parameter 5 S for a fixed value of 5 C = 5 x 10 -4 . In part b) the total 
NLO cross section is plotted versus the collinear cutoff parameter 5 C for a fixed value of 5 S = 10 -2 . 
The 4- and 5-body contributions are also shown. The cuts listed in Sec. IIIB have been imposed. 

FIG. 2. The invariant mass of the four leptons in the process pp — > ZZ + X — > e~e + /i~/x + + X. 
Parts a) and b) are for the Tevatron and LHC center of mass energies, respectively. The NLO 
(solid line) and LO (dashed line) cross sections are shown. The cuts listed in Sec. IIIB have been 
imposed. 

FIG. 3. Same as Fig. [2] but for the inclusive lepton transverse momentum distribution. The 
0-jet (dotted line) and Tjet (dot-dashed line) exclusive cross sections are also shown. 

FIG. 4. Same as Fig. [2] but for the inclusive lepton pseudorapidity distribution. 

FIG. 5. Same as Fig. [2] but for the angular distribution of the e~ . The angle 9 e - is measured in 
the parent Z-boson rest frame with respect to the parent Z-boson direction in the ZZ rest frame. 

FIG. 6. Same as Fig. || but for the angular correlation between the lepton decay planes [see 
Eq. ©]. 

FIG. 7. The transverse cluster mass distribution for the process 
pp — ► W~W + + X — > e~i> e u e e + + X. Parts a) and b) are for the Tevatron and LHC center 
of mass energies, respectively. The NLO (solid line) and LO (dashed line) cross sections are shown. 
The cuts listed in Sec. IIIB have been imposed. 

FIG. 8. Same as Fig. ^ but for the inclusive charged lepton transverse momentum distribution. 
The 0-jet (dotted line) and 1-jet (dot-dashed line) exclusive cross sections are also shown. 
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FIG. 9. Same as Fig. [?] but for the missing transverse momentum distribution. 



FIG. 10. Same as Fig. |?| but for the inclusive charged lepton pseudorapidity distribution. 

FIG. 11. The transverse cluster mass distribution for the process 
pp } — * W + Z + X — > v e e + fj," fj, + + X. Parts a) and b) are for the Tevatron and LHC center 
of mass energies, respectively. The NLO (solid line) and LO (dashed line) cross sections are shown. 
The cuts listed in Sec. IIIB have been imposed. 

FIG. 12. Same as Fig. |ll| but for the inclusive charged lepton transverse momentum distribu- 
tion. The 0-jet (dotted line) and 1-jet (dot-dashed line) exclusive cross sections are also shown. 



FIG. 13. Same as Fig. 11 but for the missing transverse momentum distribution. 



FIG. 14. Same as Fig. O but for the inclusive charged lepton pseudorapidity distribution. 



FIG. 15. Same as Fig. 11 but for the angular distribution of the e. The angle 9 e + is measured 



in the VF-boson rest frame with respect to the VF-boson direction in the WZ rest frame. 

FIG. 16. Same as Fig. |ll| but for the angular distribution of the pT . The angle 9„- is measured 
in the Z-boson rest frame with respect to the Z-boson direction in the WZ rest frame. 
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